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COMPUTING THE EULER NUMBER OF A BINARY IMAGE 

Background 

This invention relates to image processing and, more particularly, to 
computing an Euler number for an image. 
5 Image processing refers to the analysis and manipulation of images using 

a computer. The discipline encompasses a wide array of techniques, including 
pattern recognition and other image analysis, compression and encoding, such 
as for image transmission, and image construction, to name but a few examples. 
Image processing is used in diverse fields such as astronomy, medicine, military 
10 operations, communications, geology, meteorology and so on. 

A binary or digital image is typically represented by one or more 
topological properties. Topological properties represent the geometric shape of 
the image. Where changes are made to the image itself, such as stretching, 
deformation, rotation, scaling, translation, or other rubber-sheet transformations, 
15 the topological properties do not change. These properties, therefore, are quite 
useful in image characterization. Topological properties, for example, may be 
helpful in performing shape matching, object and pattern recognition, image 
data retrieval, and so on. 

One topological property of a digital image is known as an Euler number. 
20 An Euler number is defined as the difference between the number of connected 
components of the binary image and the number of holes in the image. Like the 
other topological properties, the Euler number remains invariant despite the 
above-mentioned transformation of the binary image. 

For some classes of digital images, Euler numbers have strong 
25 discriminatory power. In other words, once the Euler number for a particular 
digital image is known, the digital image may readily be distinguished from other 
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digital images in its class. This implies that the Euler number may be used for 
more efficient searching or matching of digital images. 

For example, the Euler number may be used in medical diagnosis such as 
the detection of malaria infected cells. By calculating the Euler number of each 
5 cell image, the malaria infected cells may be identified, as the Euler number of 
an infected cell is often different from that of an uninfected cell. The Euler 
number may also be used for image searching, such as in a database of logo 
images. Critical image processing applications such as these involve large 
amounts of data and, at the same time, expect quick real-time response. Fast 
10 computation of the Euler number of an image, for some applications, is, 
therefore, an indispensable task. 

Obtaining the Euler number for a digital image, however, may be 
computationally expensive. For a pixel matrix including 150,000 pixels, for 
example, over 300,000 pixel accesses and 50,000 or more processor cycles may 
15 be expected in computing the Euler number. 

Thus, there is a continuing need for better ways to compute the Euler 
number for a binary image. 

Brief Description of the Drawings 
Figure 1 is a block diagram of a system for computing the Euler number of 
20 a binary image according to one embodiment; 

Figure 2 is a diagram of an N x M pixel matrix according to one 
embodiment; 

Figure 3 is a diagram of a binary image in the N x M pixel matrix 
according to one embodiment; 
25 Figures 4A-4D are diagrams illustrating the union and intersection of 

binary images according to one embodiment; 
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Figure 5 is a flow diagram for sequentially calculating the Euler number of 
a binary image according to one embodiment; 

Figure 6 is a flow diagram illustrating parallel calculation of the Euler 
number of a binary image according to one embodiment; 
5 Figure 7 is a logic diagram of a run processor used to compute the Euler 

number according to one embodiment; 

Figure 8 is a logic diagram of a neighboring run processor used to 
calculate the Euler number according to one embodiment; 

Figure 9 is a sample circuit for performing parallel calculation of the Euler 
10 number for a 4 x 5 image according to one embodiment; and 

Figure 10 is a block diagram of a processor-based system including Euler 
number calculation capability according to one embodiment. 

Detailed Description 
In accordance with the embodiments described herein, an Euler number 
15 representation of a binary image is calculated. The calculations may be 
performed using discrete logic components to analyze the binary image. 

In Figure 1, a system 300 includes logic for computing an Euler number to 
represent a binary image, according to principles described below. In the 
following, a binary image 16 is analyzed by the system 300 and, based upon the 
20 analysis, an Euler number 90 is computed. 

The system 300 includes circuitry for calculating the Euler number 90 of 
the binary image. According to the embodiments described below, a portion of 
the binary image 16 is analyzed such that Euler number 90 for the portion, or 
partial binary image, may be computed. Additional portions of the binary image 
25 16 are likewise analyzed, in an iterative manner, until an Euler number 90 of the 
entire binary image 16 may be obtained. 
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The Euler number 90 may be very useful in characterizing a binary, or 
two-tone, image. Such characterization, for example, may be useful in 
performing database retrieval of binary images. Further, the principles embodied 
herein may be applied to gray-tone or color images which are divided into two- 
5 tone portions. Before describing the features of the system 300 further, an 
explanation of Euler numbers, as well as some properties of Euler numbers, 
follows. 

In Figure 2, according to one embodiment, a pixel matrix 20 of size NxM, 
in which N=14 and M=10, is shown. The pixel matrix 20 includes NxM pixels 
10 12 in which each pixel 12 may have a value of either a "1," for object pixels, or a 
"0," for background pixels. 

In arriving at the Euler number 90, an analysis of the binary image 16 
within the pixel matrix 20 is made. The Euler number is the difference between 
the number of connected components, or objects, and the number of holes in 
15 the binary image 16. A simple, geometric binary image 16 is depicted in Figure 
3. 

To determine the number of connected components and holes in the 
binary image 16, a pixel 12 is analyzed in terms of adjacent pixels, or its 
neighbors, of variable granularity. For a pixel 12 at location (x,y) in the pixel 
20 matrix 20, a 4-neighborhood 14 includes the set of pixels at locations [(x-l,y), 
(x+l,y), (x,y-l), and (x,y+l)]. For example, in Figure 2, a 4-neighborhood 14 of 
pixel locations [(2,2), (4,2), (3,1), and (3,3)] surrounds the pixel 12 at location 
(3,2). 

Similarly, an 8-neighborhood 18 of pixel locations [(8,5),(8,6),(8,7),(9,5), 
25 (9,7), (10,5), (10,6), (10,7)] is shown in Figure 2, for the pixel 12 at location 
(9,6). As expected, the 8-neighborhood 18 includes eight pixels 12 surrounding 
the pixel 12 that is analyzed. The 8-neighborhood 18 of a pixel 12 at location 
(x,y) includes the set of pixels at locations [(x-l,y-l), (x-l,y), (x-l,y+l), (x,y-l), 
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(x,y+l), (x+l,y-l), (x+l,y), and (x+l,y+l)]. Neighborhoods of different sizes 
may similarly be defined. 

In Figure 3, the binary image 16, a diamond-like geometric pattern with a 
hole in its center, occupies the pixel matrix 20. In the binary image 16, a 
5 connected component is a set of object pixels, e.g., pixels with a value of "1," 
such that any object pixel in the set is in the neighborhood of at least one object 
pixel of the same set. 

Assume that the shaded pixels represent a pixel value of "1." All of the 
shaded pixels shown in the binary image 16 of Figure 3 are connected 
10 components for a 8-neighborhood 18. However, for a 4-neighborhood 14, the 
shaded pixels 12 at (4,2), (10,2), (4,8), and (10,8) would not be considered 
3 connected components to the set comprising the remaining shaded pixels 12. 

il Thus, the size of the neighborhood may affect which pixels 12 are connected to 

'i other pixels 12. 

ii 15 Another feature of the binary image 16 used to calculate the Euler number 

* is a hole. A hole is set of background pixels, e.g., pixels of value "0" such that 

'* any background pixel in the set is in the neighborhood of at least one 

J background pixel of the same set. Further, the entire set of background pixels is 

13 enclosed by a connected component. 

,J 20 Thus, in Figure 3, the five non-shaded pixels 12 at the center of the binary 

image 16, or pixels at positions (6,5), (7,4), (7,5), (7,6), and (8,5), together 
make up the lone hole of the binary image 16. Because a hole is enclosed by a 
connected component, pixels 12 at position (5,2) and (10,7) are not considered 
hole pixels. It may be noted that, for an 8-neighborhood object, a hole should 
25 have a 4-neighborhood relation, and vice-versa. 

Still another feature used to calculate the Euler number, according to one 
embodiment, is a run. A run, R(i), of the i th column (or row) of the pixel matrix 
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20 is defined to be a maximal sequence of consecutive l's in the r column (or 
row). 

If the binary image 16, known as I, consists of a single row or a single 
column i, the Euler number of the image 16 is the same as the number of runs, 
5 R(i), of the image 16, as shown by the following formula: 

Property A: E(I) = R(i) 
For example, for the 0th row of the pixel matrix 20 of Figure 3, the run 
R(0) is zero, since there are no consecutive l's in the row. For the first row, the 
number of runs, R(l), equals one. In the second row, the number of runs, R(2), 
10 is three, as there are three maximal sequences of Is: a single '1/ at position 
(4,2), a sequence of three consecutive l's, (from 6,2) through (8,2), and a single 
{3 '1/ at position (10,2). 

¥ ;1 Thus, from Property A, the Euler number, E(I), of image I may be 

■'5 i=f 

= y determined, where I includes a single row (or column). Additionally, the Euler 

": sr 

|g 15 number satisfies the additive set property. Given two images, Ii and I 2 , the 
l J union (u) of the two images is defined as a simple juxtaposition of Ii and I 2 , 

! ° either vertically or horizontally, without any overlap. 

J a For example, an image, Ii, in Figure 4A includes pixel rows 1 through k. 

O Image I 2 , in Figure 4B, includes rows 1 through m. The union of these images, 

'° 20 Ii u I 2 , is illustrated in Figure 4C. Notice that the resulting image Ii u I 2 
includes k+m pixel rows. 

In contrast, the intersection (n) of Ii and I 2 is the image formed by the 
last row (or column) of Ii, and the first row (or column) of I 2 , if the images Ii 
and I 2 are joined horizontally (or vertically). The intersection image Ii n I 2 is 
25 illustrated in Figure 4D. Note that the intersection image is always two pixel 
rows (or columns) wide. 

Suppose binary images Ii and I 2 are joined horizontally, i.e., Ii lies entirely 
above I 2 , as shown in Figure 4C. The last row of Ii lies above the first row of I 2 . 
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Two runs appearing in two adjacent rows each are said to be neighboring if at 
least one pixel 12 of a run is in the neighborhood of a pixel 12 of the other run. 
Neighboring runs are defined according to the type of neighborhood, whether an 
8-neighborhood, a 4-neighborhood, or a neighborhood of some other size. 
5 From these observations, the Euler number for an image, I, which is the 

union of partial images, Ii and I 2 , with Euler numbers E(Ii) and E(I 2 ), 
respectively, may be obtained: 

Property B: E(I) = E(Ii u I 2 ) = E(Ii) + E(I 2 ) - E(Ii n I 2 ) 
As shown in Figure 4D, the image Ii n I 2 is an image containing the last 
10 row of Ii and the first row of I 2 . No holes can be present in the two-row-wide 
image. Where 0, equals the number of neighboring runs between the (i-l)* and 
i th row, an additional observation may be made: 

Property C: E(Ii-i n 10 = Q 
So, where Ik is the partial image consisting of rows 1, 2,...,(i-l) rows of 
15 the pixel matrix 20, E(Ij-i) is the Euler number of In. The Euler number of the 
image, I, consisting of only a single row, i, is R(i), from Property A. If the row, i, 
is added to In to form the u image, Ii, the intersection of these two images is 
formed by the (i-l) th and the i* rows. 

From the above principles, an Euler number, E(I), for an image, I, may be 
20 computed incrementally from the partial images, Ii, I 2 , ... In-i, In: 
E(Ii) = R(l) 

E(l 2 ) = E(Ii) + E(2) - 0 2 = R(l) + R(2) - 0 2 
E(I 3 ) = E(I 2 ) + E(3) - 0 3 = R(l) + R(2) + R(3) - (0 2 + 0 3 ) 
and, to generalize this, the formula: 
25 E(I N ) = E(I)=E(I N -i)+E(N)-0 N 

= (R(l) + R(2) + R(3) +...+ R(N)) - (0 2 + 0 3 + ... + 0 N ) 

e(In) = SRCO-io, 
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where I N =I denotes the entire image. 

Restating the formula, the Euler number of the binary image 16 is the 
difference between the sum of the number of runs for all rows (or columns) and 
the sum of the number of neighboring runs between all consecutive pairs of rows 
5 (or columns) of the pixel matrix 20. 

Looking back to Figure 1, according to one embodiment, the system 300 
includes a run processor 70, for calculating the number of runs for the rows (or 
columns) of the binary image 16. The run processor 70 is described in more 
detail in Figure 7, below. 
10 The system 300 further includes a neighboring run processors 80, in one 

embodiment, for calculating the number of neighboring runs for consecutive 

o 

O rows (or columns) of the binary image 16. The neighboring run processor 80 is 

m described in more detail in Figure 8, below. 

j ; A math processor 66 computes the sum of the number of runs and the 

iio 15 sum of the number of neighboring runs, according to one embodiment. The 

math processor 66 additionally may subtract the two sums, as specified above, in 
H computing E(I N ). The run processor 70, neighboring run processor 80, and math 

jy processor 66 may be implemented using discrete logic, specialized processing 

H hardware such as digital signal processors (DSPs), or software, to name a few 

20 possible implementations of the system 300. 

Using the formula for calculating the Euler number, E(I N ), above, the Euler 

number of the binary image 16 may be iteratively computed. The Euler number 

of the first row (or column), comprising a first or partial image, is calculated. 

The next row (or column) is then added to the partial image and then the Euler 
25 number of the union image is calculated. Thereafter, successive rows (or 

columns) are added together and, accordingly, the Euler number of each partial 

image is calculated. The process is repeated until all rows (or columns) are 

exhausted. 




In Figure 5, a flow diagram illustrates the calculation of the Euler number 
for the binary image 16, according to one embodiment. The number of runs, 
R(l), in the first row (or column) of the binary image 60 of the pixel matrix 20 is 
computed (block 102). A first Euler number, E(I), for the partial image 
5 occupying only the first row (or column) is calculated, using property A (block 
104). The variable, i, denoting the row (or column) of the pixel matrix 20, is set 
for subsequent analysis of the binary image 16 (block 106). 

Succeeding operations are performed iteratively for all i rows (or columns) 
of the binary image 16. The number of runs, R(i), in the i th row (or column) is 
10 calculated (block 108). Likewise, the number of neighboring runs, Oj, between 
the (i-l)* and the i th rows (or columns) is calculated (block 110). 
q Once the number of runs, R(i) and the number of neighboring runs, Oj, is 

1 3 known, a new Euler number, E(I), of a partial image consisting of the first and 

J= subsequent rows (or columns) may be calculated, according to the formula 

M 15 derived above (block 112). This process is repeated iteratively until all the rows 
= ' (or columns) of the pixel matrix 20 have been analyzed. 

J* In a second embodiment, multiple run and neighboring run number 

J U calculations are performed in parallel, as shown in Figure 6. Although the binary 

o image 16 is subdivided into partial images according to rows, the binary image 

20 16 may instead be divided into partial images comprising columns. 

While the number of runs in the 1 st row is calculated, the numbers of runs 
in the 2 nd row may be calculated, the number of runs in the 3 rd row may be 
calculated, and so on up to the row, where N parallel processors are available 
(block 204). 

25 Likewise, the number of neighboring runs between the 1 st and the 2 nd 

rows may be calculated, while, simultaneously with the operations of block 204, 
the number of neighboring runs between the 2 nd and 3 rd rows, between the (N- 
1)* and N 01 rows, and so on are calculated, where (N-l) neighboring run 
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processors are available (block 206). Where N-l processors are not available, a 
smaller number of parallel operations may be performed. The variable, i (shown 
in block 202), could be incremented based upon the number of parallel 
operations being performed. The operations of blocks 204 and 206 may be 
5 repeated, as needed, until both the number of runs and the number of 
neighboring runs for all rows have been calculated. 

In one embodiment, the runs for all the rows are added together (block 
208). Likewise, the number of neighboring runs for all pair of consecutive rows 
are added together (block 210). The Euler number 90 for the complete binary 

10 image 16 is calculated by subtracting the sum of the neighboring runs from the 
sum of the neighboring runs (block 212). 

In one embodiment, a very large scale integration (VLSI) architecture is 
implemented to compute the Euler number 90 for the binary image 16. A special 
purpose VLSI chip may serve as a coprocessor to a host computer, for example, 

15 to expedite computation of the Euler number 90. See, for example, Figure 10, 
below. 

In one embodiment, the VLSI architecture includes the two processing 
elements shown in Figure 1. The run processor 70 computes the number of runs 
in a row (or column). The neighboring run processor 80 computes the number 
20 of neighboring runs between two consecutive rows (or columns). 

In Figure 7, a run processor 70, according to one embodiment, includes 
discrete elements for identifying runs within the binary image 16. The discrete 
elements include a D flip-flop 34, an AND gate 44, and a modulo counter 42. 

In one embodiment, the run processor 70 calculates the number of 
25 transitions from '0' to '1' in the N x M pixel matrix 20. The number of runs is 
equal to the number of such transitions. An input line 30 is fed into the D flip- 
flop 34, for receiving each pixel in a row (or column) of the pixel matrix 20. 
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The D flip-flop 34 is set to 0 at the beginning of processing each row (or 
column). The D flip-flop 34 holds the value of the previous pixel to check for a 
transition. The pixels in a row (or column) are pipelined into the run processor 
70, at the input line 30. At any instant of time, tj, the I th and (i-l) th pixel of a 
5 particular row (or column) are checked for a 0-to-l transition. 

In this manner, the run processor 70 identifies runs for each row (or 
column) of the pixel matrix 20. The D flip-flop 34 further receives a clock 32 to 
synchronize with each input pixel. When the clock 32 is high, a Q output 36 
follows the input line 30. When the clock 32 is low, there is no change to the 
10 output. The Q output 38 AND-ed with the input line 30 produces an output 40, 
shown as X, which, if high, identifies the start of a run for the row (or column) 
being analyzed. As shown below, the Q output 36 from the D flip-flop 34, and 
the output X from the AND gate 44, may be fed into a neighboring run processor 
80. The X output 40 is fed into the modulo counter 42, which keeps a count of 
15 the number of times X goes high, i.e., the number of runs. The X output 40 may 
also be received by the neighboring run processor 80, as shown below. 

For a binary image 16 in an N x M pixel matrix 20, if the pixels are 
processed column-by-column (or row-by-row), the maximum number of runs 
possible in a row (or column) is TM/21 (or TN/21). If the processing is done 
20 column-by-column, the counter is a modulo TM/21 counter 42. The modulo 
TM/21 counter 42 also receives the clock 32, and, each time a run for the column 
is identified, the counter 42 is incremented. An output signal 46 indicates the 
number of runs counted. 

Thus, in one embodiment, after all pixels have been processed and all 
25 runs have been counted, the output signals 46 are fed to a run adder 72, as 
shown in Figure 9, for computing a sum of all such occurrences of runs. 

The neighboring run processor 80 identifies the neighboring runs between 
consecutive rows (or columns) of the binary image 16. The neighboring run 
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processor 80 checks for the condition when a particular run in a row (or column) 
starts and also determines whether the run is in the neighborhood of another run 
in its adjacent row (or column). In one embodiment, the pixels 12 of two 
adjacent rows (or columns) are fed to the neighboring run processor 80 in 
5 pipeline. 

Additionally, the neighboring run processor 80 receives data from the 
output of the D flip-flops 34 (Q outputs 36a and 36b) and AND gate 44 (X 
outputs 40a and 40b) of two consecutive rows (or columns). Accordingly, the 
neighboring run processor 80 is flanked on either side by a run processor 70, in 
10 one embodiment. At any instant of time, ti the i th and (i-l) th pixels 12 of 
consecutive rows (or columns) are checked for a neighboring run. 
j As with the run processor 70, the neighboring run processor 80, in one 

■"I embodiment, uses discrete logic to perform counting operations. AND gates 52, 

is* 

"U 54, 56, and 58 receive the X signals 40a and 40b, the Q signals 36a and 36b, 

IS 15 and the input signals 30a and 30b, as depicted in Figure 8. Output signals 92 - 

'il 

98, received from the AND gates 52 - 58, respectively, go high when neighboring 
'* runs are identified. The output signals 92 - 98 are fed into a four-input OR gate 

U 50, to produce output signal 48. Each time the output signal 48 goes high, a 

:j neighboring run has been identified in the pixel matrix 20. 

* 20 The neighboring run processor 80 also includes a modulo counter 60, for 

counting the occurrences where the output signal 48 goes high. The modulo 
counter 60 receives the clock signal 32. 

For an N x M pixel matrix 20, the maximum number of neighboring runs is 
2xTM/2l, or M+l. Accordingly, in one embodiment, the modulo counter 60 is a 
25 modulo (M+l) counter. The counter produces output signals 62, which comprise 
a sum of the neighboring runs in a pair of consecutive rows (or columns). 

To process an N x M image in parallel, according to one embodiment, N 
copies of the run processor 70 and N-l copies of the neighboring run processor 
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80 are used. For example, in Figure 9, a circuit for calculating the Euler number 
for a 4 x 5 image includes four run processors 70 and three neighboring run 
processors 80. 

Each neighboring run processor 80 receives three sets of inputs, according 
5 to one embodiment. Two input signals 30 correspond to each pixel in adjacent 
rows (or column) of the pixel matrix 20. Two X outputs 40 and two Q outputs 36 
correspond to the neighboring relation of a run for the adjacent rows (or 
columns). 

In one embodiment, the outputs 46 from the run processors 70 are fed 
10 into a run adder 72. Likewise, the outputs 62 from the neighboring run 
processors 80 are fed into a neighboring run adder 74. The run adder 72 is an 
N-input adder and the neighboring run adder 72 is a (N-l)-input adder. In one 
embodiment, the adders 72 and 74 may perform additions up to [N x Tm/21] and 
[2 x (N - 2) x Tm/21], respectively. 
15 The run adder 72 produces output signal 82, shown as X, while the 

neighboring runs adder 74 produces output signal 84, shown as Y. The circuit of 
Figure 9 further includes a subtracter 76, for subtracting the sum of the 
neighboring runs (or Y 84) from the sum of the runs (or X 82). In one 
embodiment, the subtractor 76 handles values in the range [{-2 x (N-l) x 
20 rM/2]}, {N x TM/21}]. A result 90 is the Euler number of the 4 x 5 image 
received. 

The system 300 of Figure 1 may be implemented in a processor-based 
system 400, as in Figure 10. A processor 24 is connected to the system 300 by a 
bus 26. The system 300 may serve as a coprocessor to the processor 24, to 
25 compute the Euler number 90. 

The processor-based system 400 further includes a storage 22, such as a 
memory or other volatile medium. As examples, the storage 22 may be a 
random access memory (RAM) or a flash/read-only memory (ROM), although this 
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represents but two of many possible types of storage. Additionally, the storage 
22 may store a software program 100, such as for implementing the operations 
of Figure 5. 

In one embodiment, the storage 22 may be N x M bits to store all the 
5 pixels 12 of the binary image 16 in a N x M pixel matrix 20 (since each pixel 12 
in a binary image is either 0 or 1). To calculate the number of neighboring runs 
between two consecutive rows (or columns), the column (or row) number 
identifying the starting point and the continuity of a particular run may be 
ascertained from output values Q 36 and X 40 (Figure 7). 
10 A worst case scenario for calculating the Euler number 90 arises when 

every alternating pixel 12 is an object pixel. There, the maximum number of 
3 runs in a row (column) is TM/21. Where the rows (columns) are analyzed 

J sequentially, such as depicted in Figure 5, the number of neighborhood runs is 

'i calculated between two consecutive rows (columns) at a time. In one 

8 15 embodiment, the space for storing a run is 2 locations. Under worst case 
conditions, therefore, the space for storage is 2 x 2 x TM/21. Thus, according to 
* one embodiment, the total storage space used is N x M + 4 x TM/2l. 

U The time complexity of calculating the Euler number 90 according to the 

3 embodiments described herein may be measured in terms of the number of 

20 accesses of the pixels 12 of the pixel matrix 20. In one embodiment, the 
number of pixel accesses for computing the Euler number 90 for a binary image 
16 is: 



(NxM) + 2x 



( N 



5X0 + 10, 



V i=l i=2 J 

The number of pixel accesses depends on the run distribution of the binary 
25 image 16 and, for a substantial number of binary images, R(i) << M and Oj << 
M. 
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While the present invention has been described with respect to a limited 
number of embodiments, those skilled in the art will appreciate numerous 
modifications and variations therefrom. It is intended that the appended claims 
cover all such modifications and variations as fall within the true spirit and scope 
of this present invention. 
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